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In this paper, we investigate how much of the numerical artefacts introduced by finite 
system size and choice of boundary conditions can be removed by finite size scaHng, 
for strongly-correlated systems with quasi-long-range order Starting from the exact 
C/2 . ground-state wave functions of hardcore bosons and spinless fermions with infinite 

nearest-neighbor repulsion on finite periodic chains and finite open chains, we compute 
the two-point, density-density, and pair-pair correlation functions, and fit these to vari- 
ous asymptotic power laws. Comparing the finite-periodic-chain and finite-open-chain 
rj^ . correlations with their infinite-chain counterparts, we find reasonable agreement among 

P3 ' them for the power-law amplitudes and exponents, but poor agreement for the phase 

O . shifts. More importantly, for chain lengths on the order of 100, we find our finite-open- 

chain calculation overestimates some infinite-chain exponents (as did a recent density- 
matrix renormalization-group (DMRG) calculation on finite smooth chains), whereas 
Cn ' our finite-periodic-chain calculation underestimates these exponents. We attribute this 

systematic difference to the different choice of boundary conditions. Eventually, both 
finite-chain exponents approach the infinite-chain limit: by a chain length of 1000 for 
p^ . periodic chains, and > 2000 for open chains. There is, however, a misleading apparent 

CN| ' finite size scaling convergence at shorter chain lengths, for both our finite-chain expo- 

rt — ^ . nents, as well as the finite-smooth-chain exponents. Implications of this observation 

(^ ' are discussed. 
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H ' !• Introduction 



Long-range and quasi-long-range orders in strongly-correlated systems are the ob- 
ject of numerous theoretical yj, |2|] and numerical IJ-lsl] studies. While it is fairly 
straightforward to ascertain the presence or absence of true long-range order [6U10], it 
is much harder to identify the dominant correlations from amongst several competing 
quasi-long-range orders in the quantum-mechanical ground-state wave function. Much 
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of this difficulty lies with the fact that few models that describe strongly-correlated 
systems can be solved analytically, and we have to resort to approximate or numerical 
methods to solve for the ground-state wave function. In particular, numerical solu- 
tions can only be obtained for a finite system, when our real interest is in the dominant 
quasi-long-range order of the infinite thermodynamic system. An obvious solution to 
this finite size problem would be to harness as much computational power as we can, 
to simulate extremely large systems lllll4l3il . Even so, our state-of-the-art simulations 
are still many orders of magnitude smaller than systems probed experimentally. Alter- 
natively, we can simulate systems of different sizes, and thereafter perform finite size 
scaling il4l4l6ll . to extrapolate the results to infinite system size. 



In finite size scaling, the computational condensed matter community is essentially 



guid ed by general results from statistical mechanics Ill7l - l20ll . as well as exact solutions 



I 2ll423[l . However, finite size scaling is still very much an art, with no rigorous theorem 
proving that it will always work. Besides the finite number of sites, our choice of 
boundary conditions might also affect the convergence properties of finite size scaling, 
and hence the reliability of our numerical solutions. In two recent studies, various 



quasi-long-range correlation functions were calculated using exact diagonalization 124 1 



and DMRG II25I1 . In particular, in Ref. ||25|], the power-law exponents obtained after 



finite size scaling do not agree with those derived from exact solutions in Ref. [26]. 

In this paper, we investigate the relative importance of finite system size and choice 
of boundary conditions, by comparing the power-law exponents of various correlation 
functions in a finite periodic and open chain of hardcore bosons or spinless fermions 
with infinite nearest-neighbor repulsion, against those of the infinite chain,|l26!]_as well 
as exponents of finite smooth chain obtained from the DMRG calculation [25]. In 
Sec.|2l we will describe how the ground state of finite periodic/open chain of particles 
with infinite nearest-neighbor repulsion (excluded chain) can be mapped to the ground 
state of finite periodic/open chain of particles without nearest-neighbor repulsion (in- 
cluded chain). Based on this correspondence between ground states, we then explain 
how correlation functions in the excluded chain can be written in terms of correspond- 
ing correlation functions in the included chain, through an intervening-particle expan- 
sion. Thereafter, we present in Sec. [3] our results for the two-point, density-density, 
and pair-pair correlation functions for finite excluded chains of hardcore bosons and 
spinless fermions. We fit these to different asymptotic combinations of simple and 
oscillatory power laws, to find generally good agreement between the amplitudes and 
exponents, and generally poor agreement between the phase shifts for finite and infinite 
chains. More importantly, we find for finite chains of the same intermediate lengths, 
reasonable agreement between our open-chain exponents and the finite-smooth-chain 
exponents in DMRG study. Both show strong systematic differences with the finite- 
periodic-chain exponents. By repeating these exponent calculations for very long 
chains, we established that these differences at intermediate chain lengths to be due 
to the different choice of boundary conditions. These findings are then summarized in 
SecS 



2. Models and Methods 

This section is organized into two subsections. In Sec. 12.11 we define our chain 
models of hardcore bosons and spinless fermions with infinite nearest-neighbor repul- 
sion, and show how the ground states of these finite periodic/open excluded chains 
can be mapped to the ground states of finite periodic/open included chains of hard- 
core bosons and spinless fermions without nearest-neighbor repulsion through a right- 
exclusion map. Then in Sec. 12. 21 we describe the intervening-particle expansion method, 
which allows us to write the ground-state expectation of a given excluded-chain ob- 
servable, as a sum over conditional ground-state expectations of included chain ob- 
servables. 

2.1. Chain models 

The excluded chain models for hardcore bosons and spinless fermions are given by 
the Hamiltonians 

J J 

Hie.f) ^ _f ^ (cJQ„ + Cl.Cj) + y 2 NjNj.u 
J J 

respectively, where t is the hopping matrix element, V is the nearest-neighbor repulsion, 
and U is the on-site repulsion for hardcore bosons. Exact solutions for these models 
can be obtained in the limit f/ — > oo and V -* oo, such that each site can only be singly 
occupied, and no adjacent sites can be simultaneously occupied. 

The excluded chains can be mapped to included chains of hardcore bosons and 
spinless fermions, with Hamiltonians 

//('» = -tj]{b]bj., + bX^^bi) + uY,nj(nj - 1), 

' (2) 



i/<^./).-r2(c]c,,,+c;,c,). 



respectively, using a right-exclusion map first used by Fendley to map a supersym- 
metric chain of spinless fermions to the XXZ chain 1.30.1 . As shown in Fig. [U the 
right-exclusion map deletes an empty site to the right of each occupied site, to map 
an excluded configuration for a chain of length L to an included configuration for a 
chain of length L- P. For finite open chains, this mapping is one-to-one, as shown in 
Ref. 12611 . For finite periodic chains, the right-exclusion map is not one-to-one, as illus- 
trated in Fig. 1. Nevertheless, it is possible to construct a one-to-one correspondence 
between the Bloch states of the excluded chain and Bloch states of the included chain. 
Details for this construction can be found in Ref. 13111 . 

For the rest of this paper, we will use L to denote the length of the periodic/open 
excluded chain, whose sites are indexed hy j - 1 , . . . , L. We will also use P to denote 
the total number of particles on the chain, and Aj - Bj, Cj as a common notation for 
hardcore boson and spinless fermion operators. 
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Figure 1 : Schematic diagram illustrating the many-to-one mapping from P-particle configurations of an 
excluded periodic chain of length L to P-particle configurations of an included periodic chain of length 
L' = L- P,hy deleting a single empty site on the right of a particle. 



2.2. Intervening-particle expansion 

Because of the one-to-one correspondence between Bloch states, the Bloch-state 
ampUtudes in the excluded chain ground state are identical to those in the included 
chain ground state. Hardcore bosons or spinless fermions, the included-chain ground 
state can be ultimately written in terms of a one-dimensional Fermi sea with discrete 
wave numbers. Hence the ground-state expectations of all observables on the excluded 
chain can ultimately be expressed in terms of such Fermi-sea expectations. In this 
subsection, we briefly describe how the correlations (O1O2) of two separated local 
operators of the excluded chain can be calculated using the method of intervening- 
particle expansion. 

First, for every observable O on the excluded chain, we note that it is possible to 
define a corresponding observable O' on the included chain, such that 



i(rioi7^> = i</lo'i/), 

A? n 



(3) 



for all excluded configurations \]") and \fi) that are mapped to the included configura- 
tions \f) and \f) by the right-exclusion map. Eq. Q then tells us that 



A? n 



(4) 



between the ground-state expectations of O and O' . 



Next, we expand the ground-state expectation 

{OjOj„-)^J]<^jO,Oj„.} (5) 

ip] 

as a sum over conditional ground-state expectations, where Op is a product of p particle- 
occupation number operators A^^ and r - p hole-occupation number operators (1 - Nj). 
We call this sum over all possible ways to insert particles between Oj and Oj+r the 
intervening-particle expansion. Finally, we rewrite each conditional ground-state ex- 
pectation (OjOpOj+r) on the excluded chain in terms of their corresponding condi- 
tional ground-state expectation ((9'(9' O'^^. ) on the included chain. In terms of these 
conditional expectations, which can in turn be written in terms of expectations of the 
one-dimensional Fermi sea, the intervening-particle expansion becomes 



<o,o,-..> = fX<o;o;o;.,,_p, (6) 



{p\ 
which can always be evaluated numerically for separations r that are not too large. 



3. Results and discussions 

In Ref. [26], three correlation functions: (i) the two-point function (A -Aj+r); (ii) the 
density-density correlation function (NjNj+r)', and (iii) the pair-pair correlation func- 
tion {A.A.^2^j+rAj+r+2} were systematically examined for infinite chains of hardcore 
bosons and spinless fermions with infinite nearest-neighbor repulsion, and also in three 
limiting cases for an infinite ladder. In this section, we calculate these three correlation 
functions for the finite periodic/open chain, for comparison against those of the infi- 
nite chain rt26ll . as well as against those of finite smooth chain computed in the recent 
DMRG study pS"). For the finite open chain, we analyze the numerical correlation 
functions only for N < 0.25. We believe the results are not reliable for N > 0.25, be- 
cause of phase separation in the finite open chain ground states (the symptom of which 
can be seen in Fig. |7]i. 

To understand the relative importance of finite chain length and choice of boundary 
conditions, we fit all the slowly-decaying correlation functions to the asymptotic form 

Ao -H A 1 r-"^- + A3 r-^* cos(kr + (f>), (7) 

where Aq is a constant term, Ai and A2 are the amplitude and exponent of a simple 
power law, while A3 and A4 are the amplitude and exponent of an oscillatory power 
law with wave number k and phase shift (f>. We then check how these finite-chain 
exponents depend on chain length and boundary conditions, and thus the efi'ectiveness 
of finite size scaling. 

3.1. Two-point functions 

The hardcore boson two-point function {B'.B j+i) obtained using the intervening- 
particle expansion appears to be an oscillatory power law sitting on top of a simple 
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Figure 2: Two-point function {B Bj-f-,-) of liardcore bosons witli infinite nearest-neighbor repulsion on (a) 
a periodic and (b) an open chain of length L = 100 as a function of the separation r. In this figure, red 
circles, blue aquares, and black triangles are the exact correlations at densities (a) N = 0.20, 0.30, 0.40 and 
(b) N = 0.05, 0.10, 0.15, respectively. The red, blue and black curves are the best nonlinear fits of the form 

A[ r^'^2 -I- /ij i-^'^t cos(kr + (p) to the exact correlations at these densities. 



power law. As shown in Fig. |2] we fit the hardcore boson two-point function to the 
asymptotic form {B .Bj+,} - Ai r '*- +A3 r^^* cos(kr + <f>). As expected, the fitted wave 

number is fc = ttN. 

For the infinite chain, the leading exponent (that of the simple power law) was 
found to be very close to A2 - 5, as predicted by Efetov and Larkin for an included 
chain of hardcore bosons ll29ll . For our finite periodic chain, this leading exponent 
was found to approach | with increasing density from below, as shown in Fig. |3] In 
comparison, the finite-smooth-chain exponent A2 - 2N + j obtained in Ref. [25] ap- 
proaches 5 with decreasing density from above. The finite-open-chain exponent A2 
agrees very closely with the infinite-chain results. Otherwise, we find reasonable agree- 
ment between the amplitudes and exponents of the finite and infinite chains, but poor 
agreement between their phase shifts. 

Similarly, for the spinless fermion two-point function, which is fitted to a single 
oscillatory power law {C -Cj+r) - A^ r"'*-' cosikr + (p) (see Fig.|4|, we find reasonable 
agreement between the amplitudes and exponents of the finite and infinite chains, but 
poor agreement between the phase shifts of the finite and infinite chains (see Fig. |5]l. 
As expected, the fitted wave number is A: = nN. For noninteracting spinless fermions, 
the exponent of the oscillatory power law is A4 = 1. We see from Fig.|5]that for both 
finite and infinite chains, A4 starts close to 1 at A^ = 0, and decreases with increas- 
ing density. In particular, the finite-chain exponent is smaller than the infinite-chain 
exponent, telling us that two-point correlations are stronger on finite chains. 



3.2. Density-density correlations 

The density-density correlation function, which is identical for hardcore bosons and 
spinless fermions on a chain, can be fitted very well to the asymptotic form {NjNj + r) = 
Aq + A3 r^'^^ cos(kr + 0) (see Fig.|6l). Again, the fitted wave number k = InN agrees 
with our expectations. Also, as with the two-point functions, we see in Fig.|2]that there 
is reasonable agreement between the finite and infinite chains for the fitted amplitudes 
and exponents, but poor agreement for the phase shifts. For both finite chains, we find 




Figure 3: The fitted (top) amplitudes A[ (black) and A3 (red), (middle) exponents A2 (black) and A4 (red), 
and (bottom) phase shift (p of the simple power law and oscillatory power law as a function of density N, 
for the two-point function of a L = 100 periodic(dashed lines) and open(dot-dashed lines) chain of hardcore 
bosons with infinite nearest-neighbor repulsion. The corresponding parameters of the infinite chain from 
Ref. I2al are shown as solid lines, whereas finite-smooth-chain results from Ref. 12511 are shown as dotted 
lines. 
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Figure 4: Two-point function {C Cj+,-) of spinless fermions with infinite nearest-neighbor repulsion on (a) 
a periodic and (b) an open chain of length L = 100 as a function of the separation r. In this figure, red 
circles, blue squares, and black triangles are the exact correlations at densities (a) N = 0.2,0.3,0.4 (b) 
N = 0.05,0.10,0.15, respectively. The red, blue and black curves are the best nonlinear fits of the form 
A3 r^^* cos(kr + (p) to the exact correlations at these densities. 
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Figure 5: The fitted amplitude A3 (top), exponent A4 (middle) and phase shift ^ (bottom) of the oscillatory 
power-law fit to the two-point function of a L = 100 periodic(dashed lines) and open(dot-dashed lines) chain 
of spinless fermions with infinite nearest-neighbor repulsion, as a function of density A'. The con'esponding 
parameters of the infinite chain from Ref. |26] are shown as solid lines. 



Aq > (Nj) {Nj + r) - N", the limit attained in the infinite excluded chain. This ex- 
cess density-density correlation compared to the infinite chain persists at L = 200 and 
L - 1000. Also shown in Fig.|2]are the oscillatory power- law exponent A4 for the finite 
periodic/open chain, the infinite chain, and the finite smooth chain. 

In Fig. 13 a), while we see the exponents in all three cases decaying with increas- 
ing N, we again find the finite-periodic-chain exponent (dashed line) to be generally 
smaller than the infinite-chain exponent (solid line), the finite-open-chain exponent to 
be slightly larger than the infinite-chain exponent, and the finite-smooth-chain exponent 
(dotted line) to be larger than the infinite-chain exponent. Since the same chain lengths 
were used, we suggest two possible explanations for these systematic differences. As 
we shall later show for these intermediate chain lengths, the weaker density-density 
correlation of finite smooth chain and finite open chain is not an artefact resulting from 
spatial averaging done to 'restore' translational symmetry. Instead, it is a genuine fi- 
nite size effect related to the choice of smooth [25] and open boundary conditions, as 
opposed to periodic boundary conditions used on our finite periodic chains. 



3.3. Pair-pair correlations 

The pair-pair correlation function (A . jA -Aj+rAj+r+i) is also identical for both spin- 
less fermions and hardcore bosons. For spinless fermions, {A ._2^Aj+rAj+r+2} can be 
interpreted as the superconducting correlations. As shown in Fig.[8J the pair-pair corre- 
lation function can be fitted very well to the asymptotic form Air"'*^ +Air^'^-' cos{kr+(p). 
As expected, the wave number isk = 2nN if we allow A; to be a fitting parameter. Again, 
from Fig.|9l we see that there is reasonable agreement on the amplitudes and exponents 
between the finite and infinite chains, but poor agreement between their phase shifts. 
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Figure 6: Density-density coirelation function (NjNj+r) of particles witli infinite neai'est-neiglibor repulsion 
on a periodic/open chain of length L = 100 as a function of the separation r. In this figure, red circles, 
blue squares, and black triangles are the exact correlations at densities (a) N = 0.15,0.30,0.45, and (b) 
N = 0.05,0.10,0.15, respectively. The red, blue and black curves are the best nonlinear fits of the form 
Aq + A3 r"''-' cos(A:r + <f>) to the exact correlations at these densities. 




Figure 7: The fitted constant Aq (top), and amplitude A3 (second from top), exponent A4 (second from 
bottom), and phase shift </> of the oscillatory power law as a function of density N, for the density-density 
congelation function of a L = 100 periodic(dashed lines) and open(dot-dashed lines) chain of particles with 
infinite nearest-neighbor repulsion. The cortesponding parameters of the infinite chain from Ref. 1 26] are 
shown as solid lines (A4 = ^ + j{^ - N) for the infinite chain), while the finite-smooth-chain exponent 
A4 = -2.96 N + 2.23 is shown as the dotted fine. 





Figure 8: Pair-pair correlation function (A .^A-J^j+r^j+r+l) for particles with infinite nearest-neighbor re- 
pulsion on (a) a periodic (b) an open chain of length L= 100 as a function of the separation r. In this figure, 
red circles, blue squares, and black triangles are the exact correlations at densities (a) N = 0.20, 0.25, 0.30 
(b) N = 0.05,0.10,0.15, respectively. The red, blue and black curves are the best nonlinear fits of the form 
Ai r"'*2 -I- A3 r^'^* cos{kr + ^) to the exact correlations at these densities. 




Figure 9: The fitted amplitudes Ai and A3 (top), exponents A2 and A4 (middle), and phase shifts ^ (bottom) 
of the pair-pair correlation function of a L = 100 periodic(dashed lines) and open(dot-dashed lines) chain 
of particles with infinite nearest-neighbor repulsion, as a function of the density N. The corresponding 
parameters of the infinite chain from Ref. [26] are shown as solid lines. 
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3.4. Comparison between different chain lengths 

In the DMRG study of finite smooth chain [25], a systematic study on the effects 
of chain lengths was undertaken, even though there was no rigorous attempts at per- 
forming finite size scaHng. They found, for the oscillatory power-law exponent A4 
of the density-density correlation on ladders of lengths L = 100, 150,200, the fitted 
parameter appears to have converged onto the straight line A4 = -2.96 N + 2.23. In 
Fig. [TOl we show the fitted exponent of the density-density correlations on finite pe- 
riodic chains with lengths L = 50, 100,200, 1000 and finite open chains with lengths 
L = 100, 150, 200, 300, 1000, 2000. For the finite periodic chain, we also find an appar- 
ent convergence at a chain length of approximately L = 200. However, this 'limiting' 
behaviour of A4 does not agree with the exponent determined from the infinite chain, 
and is also systematically different from the 'limiting' finite-smooth-chain exponent. 
In fact, unlike the infinite-chain and finite-smooth-chain exponents, the finite-periodic- 
chain exponent appears to be a strongly nonlinear function of N. When we go to a 
L - 1000 periodic chain, the exponent A4 is still strongly ./V-dependent. However, we 
can see from Fig. [10] that it is closer to the infinite-chain limit. 

For finite open chains, however, the chain-length-dependence of the exponent A4 
is rather different. Starting from being close to the infinite-chain limit at L = 100, we 
see the exponent moving away from the infinite-chain limit as the chain length goes to 
L = 150 to L = 200 to L = 300 to L = 1000. In fact, there is an apparent convergence 
at L = 200 to the finite-smooth-chain exponent. This suggests that a L - 100 chain 
with smoothed open boundary conditions behaves like a longer chain with hard open 
boundary conditions, which is what they desired when White introduced the smoothed 



boundary conditions II27L 12811 . What they, and DMRG group would not have guessed, 
is the finite open chain exponent deviating from the infinite chain exponent initially 
as we go to longer and longer chains. From the open chain results of L = 1000 and 
L - 2000, we find that the convergence to infinite-chain results is slower than for finite 
periodic chains. 

With this, we demonstrated that the choice of periodic/open boundary conditions 
resulted in the finite-chain exponents being systematically lower/higher than the infinite- 
chain exponent, for chain lengths on the order of L = 100 to L = 200. This observation 
is important, because DMRG calculations on finite smooth chains are currently lim- 
ited to such chain lengths. Furthermore, we showed that the finite-chain exponents do 
converge onto the infinite-chain limit, but only when the chain length is on the order 
of L = 1000. Moreover, this convergence is not uniform for finite chains. This poses 
an important challenge for the use of finite size scaling to determine the infinite-chain 
limits. 

Finally, we compare the density-density correlation functions of different chain 
lengths, against the infinite-chain density-density correlation function, at density N - 
0.10. As shown in Fig.[TT] oscillatory power-law decay in the density-density correla- 
tion function is seen for finite and infinite chains. The finite periodic chain and infinite 
chain amplitudes are very similiar. The rates at which the density-density correlations 
decay on the finite periodic chain and infinite chain are also very similar. This explains 
the reasonable agreement seen in the fitted amplitudes and exponents. 

However, the phase shifts in the finite periodic chain and the infinite chain are rather 
different, and increasing the chain length from L = 50 to L = 100 to L = 200 does not 
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Figure 10: The fitted oscillatory power-law exponent A4 for the density-density congelation function of hard- 
core particles with infinite nearest-neighbor repulsion as a function of density N, for dilferent periodic chain 
lengths, L = 50 (red), L = 100 (blue), and L = 200 (black). The dashed line is the infinite-chain exponent 
from Ref. 1231 . whereas the dotted line is the best-fit straight line to the finite-smooth-chain exponent from 
Ref. m. 



make the finite-chain phase shift approach the infinite-chain phase shift. This is true 
even at L = 1000. As pointed out by the referee, this is expected as the phase shift is a 
parameter that is adjusted globally to minimize the effects of the boundaries, and thus 
will not converge smoothly with increasing chain length. 

More importantly, the amplitude of the oscillatory power law is very small for finite 
open chains, because of the spatial averaging performed. This same spatial averaging 
is done in the finite-smooth-chain calculations. 
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Figure 1 1 : The density-density correlation functions of infinite, finite periodic, finite open chains with N - 
0.10. The solid line is the infinite-chain exponent from Ref. 12611 . 
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4. Conclusions 

In conclusion, we employed the exact mapping between the ground states of the 
excluded chains and included chains of hardcore bosons and spinless fermions, along 
with the intervening -particle expansion method, both developed in Ref. |26], to in- 
vestigate how the finite chain length and the choice of the boundary conditions man- 
ifest themselves in the two-point, density-density, and pairing correlation functions. 
We fitted these finite-chain correlations to the generic asymptotic form Aq + A[ r"'*- -H 
A3 r"^^ cos(A:r + (p), and compared the fitted parameters to those obtained from an infi- 
nite chain, and a finite smooth chain. In general, we find reasonable agreement between 
fitted amplitudes and exponents for the finite and infinite chains, but poor agreement 
for the fitted phases. We also find the phase shifts depend on both boundary conditions 
and system length. 

Comparing the finite-periodic-chain, finite-open-chain, finite-smooth-chain and infinite- 
chain exponents for correlation functions at the same length, we find that the finite- 
periodic-chain exponents are closer to the infinite-chain results and finite-open-chain 
exponents are closer to the finite-smooth-chain ones. We attribute this systematic dif- 
ference to the finite size effects and the different choices of boundary conditions. In 
a nutshell, finite-periodic-chain correlations are stronger than they are with open or 
smooth boundary conditions (or any variants thereof). More importantly, we observed 
at intermediate chain lengths an apparent 'convergence' of our finite-chain exponents 
to functions of the density A^ different from the infinite-chain exponents. A straightfor- 
ward finite size scaling may thus lead us to the wrong physics. This same phenomenon 
was observed in the DMRG study of finite smooth chain as well, suggesting in general 
that uniform convergence cannot be expected for chain lengths of L ~ 10^, and perhaps 
not even for chain lengths of L ~ 10^. A more positive take on these general results 
would that that numerical studies using more than one set of boundary conditions offer 
us a better sense of how far we might be away from the infinite-system limits, even 
when there is no guarantee that the results from different boundary conditions will 
bound the true infinite-chain limit. 
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